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ABSTRACT: Rate coefficients of the title reaction R 31 (SO 2 + O + M -> SO 3 + M) and R 56 (SO 2 + 
HO 2 SO 3 + OH), important in the conversion ofS(IV) to S(VI), were obtained at T =970-1 150 
K and p ave = 16.2 pmol cm -3 behind reflected shockwaves by a perturbation method. Shock- 
heated H 2 / 02 /Ar mixtures were perturbed by adding small amounts of SO 2 (1%, 2 %, and 3%) and 
the OH temporal profiles were then measured using laser absorption spectroscopy. Reaction 
rate coefficients were elucidated by matching the characteristic reaction times acquired from 
the individual experimental absorption profiles via simultaneous optimization of &31 and k 56 
values in the reaction modeling (for satisfactory matches to the observed characteristic times, 
it was necessary to take into account R 56 ). In the experimental conditions of this study, R 31 is in 
the low-pressure limit. The rate coefficient expressions fitted using the combined data of this 
study and the previous experimental results are ^3 1 ,o/( Ar] = 2.9 x 10 35 T _6, ° exp(— 4780 K/T ) + 
6.1 x 10 24 T _3, ° exp( — 1980 K/T) cm 6 moH 2 s ” 1 at T = 300-2500 K; k 56 = 1.36 x 10 11 
exp(— 3420 K/T ) cm 3 mol -1 s -1 at T = 970-1 150 K. Computer simulations of typical aircraft 
engine environments, using the reaction mechanism with the above 631,0 and k 56 expressions, 
gave the maximum S(IV) to S(VI) conversion yield of ca. 3.5% and 2.5% for the constant density 
and constant pressure flow condition, respectively. Moreover, maximum conversions occur at 
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rather higher temperatures (M200 K) than that where the maximum 631,0 value is located 
(~800 K). This is because the conversion yield is dependent upon not only the 631,0 and 656 
values (production flux) but also the availability of H, O, and HO 2 in the system (consumption 
flux). © 2010 Wiley Periodicals, Inc.* Int J Chem Kinet 42: 168-180, 2010 


INTRODUCTION 

Contrails and contrail-induced cirrus clouds may af- 
fect the global climate via radiative forcing and the 
alteration of atmospheric chemistry [1-3]. The critical 
chemical species in contrails, in addition to soot, ions, 
and water vapor, are sulfates and sulfuric acid formed 
by oxidation of fuel sulfur, i.e., low oxidation state 
fuel bound S(II), mostly thiols and sulfides, are rapidly 
and nearly quantitatively converted to S(IV), i.e., S0 2 , 
in the high-temperature oxidative environment of the 
combustor, and then on to S(VI), mostly SO3 and 
small amounts of sulfates, formed in the cooler, but 
still harsh environment of the turbine and exhaust noz- 
zle [4-6] . Sulfates and sulfuric acid, formed in the air- 
craft exhaust plume, enhance the formation of contrails 
and cirrus clouds through homogeneous binary nucle- 
ation of H 2 SC>4/H 2 0 and heterogeneous condensation 
of H 2 SC>4/H 2 0 on soot particles or ions [4-9]. 

Owing to the significant role of S(VI) (SO3 and 
H 2 S04) in the formation of contrail and cirrus cloud, 
the conversion fraction of S(IV) to S(VI) has been mea- 
sured or calculated by various authors: The measured 
fractions were about 0.34%-4.5% [10,11] or 3.3% ± 
1.8% [12], while calculations from measured particu- 
late volume, presuming only H 2 S04 and H 2 0 are in 
the particles, span a wide range from 0.4% to 80% 
[10,11,13-15]. 

The route to sulfuric acid formation in the homoge- 
neous gas phase has been considered as Ro (fuel-bound 
S S0 2 ), R31 (S0 2 + O + M ^ S0 3 + M), R 32 
(S0 2 + OH (+M) -> HOS0 2 (+M)), R 52 (HOS0 2 + 
0 2 -* S0 3 + H0 2 ), and R 57 (S0 3 + H 2 0 -> H 2 S0 4 ). 
It is believed that Ro occurs inside the combustor, R31, 
R3 2 , and R5 2 in the turbine and engine exit nozzle, 
and R57 in the exhaust plume. R31 is considered as the 
principal reaction pathway for the conversion of sulfur 
from the S(IV) to the S(VI) oxidation state. As seen in 
Fig. 6, the experimentally measured rate coefficients of 
the spin-forbidden R 31 (S0 2 (X 1 A 1 ) + 0( 3 P) + M -> 
S0 3 (X 1 A / 1 ) + M) exhibit an unusual behavior — an in- 
crease with temperature ( T <~450 K) [16,17] followed 
by a sharp decrease with further increase of tempera- 
ture ( T >~1400 K) [18-20]. This unusual behavior of 
rate coefficients with respect to temperature was first 
rationalized by Astholz et al. [19] in terms of potential 
energy barrier arising from the triplet-singlet intersys- 
tem crossing and diminishing collision efficiency with 


temperature in the stabilization of the collision com- 
plex formed (see the section Discussion for details). 

Despite its important role in contrail formation, the 
rates of this reaction have never been measured at tem- 
peratures typical of aircraft operation (900-1200 K). 
Moreover, there are large discrepancies among the ex- 
trapolated values from the existing rate coefficient ex- 
pressions to this temperature region. For example, at 
T = 1000 K the calculated or extrapolated values from 
the early studies are (in cm 6 mol -2 s -1 ) 2.8 x 10 15 

[21] , 4.1 x 10 15 [17], 9.6 x 10 15 [16], 1.0 x 10 16 
[20], and 3.2 x 10 16 [19]. An order of magnitude dis- 
crepancy is easily recognized. Recently, Naidoo et al. 

[22] measured the rate coefficients of R31 in the tem- 
perature range of 290-840 K and reported the Rice- 
Ramsperger-Kassel-Marcus (RRKM) fitted parame- 
ters for &3i ; o/[Ar], £31,00 and broadening factor F c for 
T = 220-2500 K. The £ 3 i, 0 /[Ar] value of Naidoo et al. 
at T = 1000 K gives 2.8 x 10 15 cm 6 mol -2 s _1 . We 
have previously reported results for this reaction over 
a more limited temperature range [23]. In the present 
work, we extended our previous efforts and completed 
the measurement of the rate coefficients of R31 at T = 
970-1 150 K using a shock tube and OH laser absorp- 
tion technique. 

SELECTION OF EXPERIMENTAL 
CONDITIONS AND COMPUTER 
SIMULATION 

Determination of the rate coefficients of R31 was per- 
formed via a perturbation method utilizing the abun- 
dant and well-characterized O-atoms provided by the 
shock-induced ignition of H 2 /0 2 /Ar mixtures. This 
permitted the use of experimentally convenient mix- 
ture dilutions and reaction time scales as well as sensi- 
tive OH absorption diagnostics in the easily accessible 
near ultraviolet wavelength regime. We perturbed the 
H 2 /0 2 /Ar reaction system with small amounts of S0 2 
and recorded the differences in temporal behavior of 
the OH absorption. In this regard, construction of a ro- 
bust reaction mechanism of the S0 2 /H 2 /0 2 /Ar system 
is essential. Our reaction mechanism is given in Table 2 
in the Supporting Information. 

For the temperature and pressure range relevant 
to S(IV) to S(VI) conversion, the major O-atom 
production reaction is Ri (H + 0 2 -> OH + O) and the 
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Sensitivity 


Figure 1 Sensitivity spectra for 150 of the S 02 /H 2 / 02 /Ar = 
2.0/0.5/10.0/87.5 mixture at T = 1036 K and p = 15.8 pmol 
cm -3 (Fig. 2 condition). 


main O-atom consumption reaction is R 2 (O + H 2 -> 
OH + H). O-atom production rates are further modified 
by R 9 (H + O 2 + M HO 2 + M) through the chan- 
neling of some H-atoms away from chain branching. In 
the S 02 /H 2 / 0 2 /Ar system, R 31 (S0 2 + O + M -> S0 3 
+ M) and R 2 compete for O-atoms. In rich H 2 / 02 /Ar 
mixtures, this competition would be dominated by R 2 
and thus the addition of SO 2 would barely affect the 
course of reaction. Therefore, we restrict our considera- 
tion to lean H 2 / 02 /Ar mixtures only. On the basis of our 
previous study [24], we selected H 2 / 02 /Ar = 0.5/10.0/ 
89.5 (0 = 0.025) as the baseline mixture. Using this 
mixture, we performed a series of sensitivity analyses 
to determine both the optimum levels of SO 2 addi- 
tion and the best set of reaction characteristics to use. 
We determined that matching the characteristic times 
(T max , T 75 , T 50 , and T 25 defined in the Results section) 
and employing three mixtures containing 1 . 0 %, 2 . 0 %, 
and 3.0% SO 2 would give the best balance of over- 
all reaction sensitivity and experimental convenience. 
Figure 1 shows sensitivity spectra obtained using the 
Table 2 reaction mechanism at a typical experimental 
condition. The sensitivity of R 31 was found to be in- 
creased with both decreasing T and increasing SO 2 /H 2 
ratio in the mixture. Although excellent R 31 sensitivi- 
ties were computed at low temperatures, the character- 
istic times became prohibitively long. 

The Table 2 mechanism was composed of 57 ele- 
mentary reactions in 18 chemical species. For H 2 /O 2 
reactions, we used our reaction mechanism of 20 ele- 
mentary reactions [24] with a change in ku (HO 2 + 


OH — H 2 0 + 0 2 ) from 2.0 x 10 13 to 3.0 x 10 13 cm 3 
mol -1 s -1 following the recommendation of Srinivasan 
et al. [25]. The H 2 /O 2 reactions were combined with 
additional 35 reactions that are related to SO 2 oxida- 
tion by H x S y O z species from Rasmussen et al. [26]. 
Then we consider the oxidation of SO 2 by HO 2 radi- 
cal, important in our experimental conditions, through 
the reaction R 56 (SO 2 + HO 2 SO 3 + OH) of which 
rate coefficients are obtained as part of the simultane- 
ous optimization with R 31 (see below). Also, we added 
R 57 (SO 3 + H 2 0 H 2 SO 4 ) [27] for final product for- 
mation for completeness. However, the contribution of 
this reaction to the reaction progress is negligible. 

All simulations were performed using our custom 
code that has advantageous built-in features for opti- 
mization. This code has been thoroughly tested against 
existing kinetic codes [28]. Recently updated values 
for the thermodynamic data [29,30] have been used 
for all species except for HOSO, HSO 2 , and HOSO 2 . 
The thermodynamic properties of these species were 
estimated using the available experimental data and 
theoretical calculations. The NASA polynomial coef- 
ficients of these species are given in Table 1 in the 
Supporting Information. 

EXPERIMENTAL 

The shock tube-laser absorption technique employed in 
this study is described in detail elsewhere [31]. Briefly, 
a rolled square stainless steel shock tube, 63.5 mm in 
cross section, was used for experiments. The test sec- 
tion was routinely pumped below 3 pTorr by a Varian 
V60 Turbopump. A combined leak and outgassing rate 
of the test section was 5 pTorr min -1 . To minimize 
problems due to wall loss of SO 2 , all fills were done 
in a two-step process. First, the shock tube was filled 
with mixture to the desired starting pressure and the 
room temperature absorption of SO 2 was measured. 
The tube was then evacuated and refilled to the same 
pressure with fresh test gas. All shocks were then ini- 
tiated within 1 min of admitting the test gas mixture. 

Shock speeds were measured using four flush- 
mounted 113A21 PCB Piezotronics pressure trans- 
ducers coupled to Phillips PM 6666 programmable 
countertimers. The last pressure transducer was axially 
coincident with the center of the observation windows, 
12.7 mm distant from the end wall. Measured shock 
speeds were fitted to a second-order polynomial in dis- 
tance. The extrapolated speed to the end wall was used 
together with NASA thermodynamic data [29] in solv- 
ing the relationships for the reflected shock properties 
[32], assuming full vibrational relaxation and no chem- 
ical reaction at the shock front. The postreflected shock 
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properties were corrected for boundary layer/reflected 
shock interaction effects using the methodology de- 
scribed by Hwang et al. [28]. 

The combined OH and S0 2 absorption profiles were 
measured at 310.032 nm (in air), corresponding to the 
P x { 5) line of the (0,0) band of the OH A 2 E + X 2 n 
transition. UV light was produced by a CW argon ion 
laser (Coherent Innova 200) pumped ring-dye laser 
(Coherent 899-21) running Kiton Red 620 dye. Intra- 
cavity frequency doubling was achieved with an an- 
gle tuned LiI03 crystal. For absorption measurements, 
a triple beam scheme was employed: a probe beam 
passed through the center of two 25.4-mm-diameter 
Sl-UVA quartz windows (12.7 mm distant from the 
shock tube end wall), an intensity reference beam, and 
a wavelength reference beam. Passage through an un- 
constrained burner- stabilized methane/air flame was 
used to provide the wavelength reference, i.e., the cen- 
ter absorption of the OH P\(5) line. Each detector was 
a combination of interference filter (X 0 = 310 ± 10 nm 
fwhm), a PMT (Thom EMI9924QB), and a high-speed 
buffer/amplifier. A custom dynode chain was used for 
these experiments to assure the linearity of the absorp- 
tion signal. The probe, reference, difference (probe de- 
reference ac), and the last pressure transducer signals 
were recorded on a four-channel Nicolet 4094C digital 
oscilloscope equipped with Nicolet 4570 plug-in units. 
Reaction progress was followed using the difference 
signal. The S/N of the difference signal was usually 
better than 60. The overall electronic time constant de- 
termined for the entire PMT-buffer-cable system was 
less than 0.5 qs. 

S0 2 /H 2 /0 2 /Ar test gas mixtures were prepared 
manometrically and allowed to stand for at least 
48 h before use. A pump-fill-pump-fill cycle was 
used for the S0 2 addition to minimize the effect 
of wall loss. The maximum uncertainty in mixture 
composition was less than 0.5% of the nominal mole 
fraction for each component. The stated purities of the 
gases were H 2 , 99.9995% (MG Industries, Malvern, 
PA; scientific grade); 0 2 , 99.998% (MG Industries; 
scientific grade); S0 2 (anhydrous), 99.98 wt% (MG 
Industries; scientific grade); Ar, 99.9999% (MG 
Industries; sputtering grade). All gases were used 
without further purification. 

RESULTS 

A typical absorption profile is shown in Fig. 2. The re- 
flected shock passage at time zero is shown in the form 
of a schlieren peak at the extreme left-hand edge of the 
experimental trace. The initial absorption is entirely 
due to S0 2 . 



Figure 2 Typical experimental and simulated absorption 
profiles. S0 2 /H 2 /0 2 /Ar = 2.0/0.5/10.0/87.5, T = 1036 K, 
and p = 15.8 pmol cm -3 . 


Because both S0 2 and OH absorb 310-nm UV 
light, the absorption coefficients of S0 2 were measured 
and used in the simulation of experimental reaction 
progress. OH absorption coefficients were calculated 
using our recent values of the spectral line broadening 
parameters and other molecular properties [33]. The 
mean room temperature base e absorption coefficient 
of S0 2 measured by pump-fill-pump-fill method was 
£ e (S0 2 ) = (7.0 it 0.4) x 10 4 cm 2 mol -1 , which is about 
30% smaller than the spectroscopic measurement of 
Vandaele et al. [34]. 

Figure 3 shows the S0 2 absorption coefficients mea- 
sured at room temperature and at temperatures behind 
incident and reflected shock waves. A temperature- 
dependent fitting gives 

s e ( SO 2 ) = 1.3 x 10 5 exp(— 180 K/7’)cm 2 m 0 r 1 . 

For reaction progress and sensitivity analyses, the fol- 
lowing set of observables was obtained from an experi- 
mental absorption profile: an initial S0 2 absorption af- 
ter reflected shock passage, Aso 2 ,o = 1 — 7so 2 ,o/G; a 
maximum absorption, A max = (1 — /// o)max; and four 
characteristic times (r max , 175, T 50 , and T 25 ). Evalua- 
tion of T max , time to reach A max , is straightforward. 
For t 2 5, however, it was necessary to use a sub- 
tracted absorption signal (A' max = A max — A S o 2 ,o) for 
the experiments that resulted in large Aso 2 ,o and small 
A max . Thus t 2 5, T50, and 175 were defined as the time 
to reach 25%, 50%, and 75% of A ' max , respectively. 
Figure 4 shows 7:50 values measured in 1 .0%, 2.0%, and 
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iooo K/r 

Figure 3 Base e absorption coefficients of SO2 at 3 10 nm. 
The least-squares fit is £ e (S02) = 1.3 x 10 5 exp(— 180K/T) 
cm 2 mol -1 . The vertical error bars represent error limits at 
the 95% confidence level. 



1000 Et /T 


Figure 4 Measured 150 values of the baseline and 
S02-perturbed mixtures. Symbols are as follows: filled 
square, S02/H2/02/Ar = 0.0/0.5/10.0/89.5 (base line mix- 
ture); star, S02/H2/02/Ar = 1.0/0.5/10.0/88.5; open cir- 
cle, S02/H2/02/Ar = 2.0/0.5/10.0/87.5; open diamond, 
S0 2 /H 2 /0 2 /Ar = 3.0/0.5/10.0/86.5. 


3.0% SO 2 /0.5% H 2 /10.0% 0 2 /Ar mixtures together 
with those from a baseline mixture (0.5% H 2 / 10 . 0 % 
02/89.5% Ar). Clearly perturbation with SO 2 prolongs 
the characteristic times. Furthermore, the perturbation 
effect increases with SO 2 /H 2 ratio and decreases with 
T . The increase of characteristic times is due mainly 


to the depletion of O atoms by R 31 , so the H-atom pro- 
duction rate by R 2 is reduced. This, in turn, reduces the 
branching rate of Ri . 

The presence of large amount of HO 2 radicals in our 
experimental systems (low temperatures, T = 970- 
1150 K and SO 2 -perturbed lean H 2 /02/Ar mixture 
compositions), convinced us to include reactions be- 
tween SO 2 and HO 2 . Three possible reactions are R _47 
(HOSO + 0 2 <- S0 2 + H0 2 ), R — 50 (HS0 2 + 0 2 <- 
S0 2 + H0 2 ), and R 56 (S0 2 + H0 2 -> S0 3 + OH). 
The first two reactions (R-47 and R- 50 ) have already 
included in the reaction mechanism of Rasmussen et al. 
[26] that was adopted in this work without any mod- 
ification. R 56 has been studied experimentally at low 
temperatures (T < 300 K) for atmospheric applica- 
tions [35-38] and theoretically at high temperatures 
for combustion applications [39,40] (see the section 
Discussion). The inclusion of the reaction R 56 , which 
has some sensitivity on the characteristic times (Fig. 1), 
is necessary to unravel the contributions by R 31 and R 56 
to the disproportionate enhancement of the character- 
istic times with [SC^Jo- 

Using the Table 2 mechanism, computer simula- 
tions for experimental observables have been per- 
formed to simultaneously determine the rate coeffi- 
cients of R 31 (SO 2 + O + M SO 3 + M) and R 56 
(SO 2 + HO 2 SO 3 + OH). Determination of the rate 
coefficients was accomplished using a method similar 
to the solution mapping [41,42], i.e., for each experi- 
ment, only the k 3 \ and k 56 were systematically varied 
and the corresponding values of the function, V(x) = 

[(^calculated (**0 ^ 1 , experimental)/^!, experimental] » where 

i = {r 2 5 , ^ 50 , T75, T max } and x = {* 31 , k 56 ], were cal- 
culated. Then an exhaustive search was made for the 
minimum value of T(x) to find the optimized values of 
x = {k 3 i , & 56 }. A simulated absorption profile is shown 
in Fig. 2. The resulting rate coefficients of R 56 and R 31 
are plotted in Figs. 5 and 6 , respectively. Experimental 
conditions, measured observables, and the rate coeffi- 
cients determined for individual experiments are given 
in Table 3 in the Supporting Information. 

DISCUSSION 

No measurements had been made for the reactions be- 
tween SO 2 and HO 2 at combustion temperatures. Three 
room temperature measurements have been made and 
calculations have been performed for combustion tem- 
peratures. Three product channels have been consid- 
ered: R _47 (HOSO + 0 2 S0 2 + H0 2 ), R — 50 (HS0 2 

+ 0 2 S0 2 + H0 2 ), and R 56 (S0 2 + H0 2 -> S0 3 + 

OH). As for the first two reactions, included in the reac- 
tion mechanism of Rasmussen et al. [26] and also used 
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Figure 5 Rate coefficients of R 56 (S0 2 ± H0 2 -> SO 3 
± OH). Symbols are as follows: circle, our data; diamond, 
Payne et al. [36]; square, Graham et al. [37]; triangle, 
Burrows et al. [38]. Lines are solid, a least-squares fit to 
our data, £55 = 1.36 x 10 11 exp(— 3420 K JT) cm 3 mol -1 
s _ 1 at T = 970-1150 K with estimated error of ±30%; 
dashed, extrapolation of our expression to room tempera- 
ture; light solid, a three-parameter fit of our data with the 
lower limit value of the room temperature data of Payne et al., 
k 56 = 1.5 x 10 " 3 7 " 4 exp(l 100 K/T) cm 3 mol -1 s ” 1 at T = 
300-2000 K; dot-dashed, calculated from /:_56 = 4.8 x 10 4 
r 2,46 exp(— 13720 K IT) cm 3 mol -1 s _1 of Hindiyarti et al. 
[39] with K c 56 . [Color figure can be viewed in the online 
issue, which is available at www.interscience.wiley.com.] 


in this study, there have been no experimental studies. 
For the third reaction channel, only the room tempera- 
ture rate coefficients are available [35-38], which vary 
from & 56 < 6.0 x 10 5 to £56 = 5.2 x 10 8 cm 3 mol -1 
s -1 . Of these, the two studies (Graham et al. [37] and 
Burrows et al. [38]) using direct in situ spectroscopic 
measurements of the reactants observed no measurable 
reaction with S0 2 and thus returned only upper limit 
values. The lower of these, Graham et al. [37], perforce 
obviates the other. Reaction with S0 2 was observed 
using the somewhat indirect photochemical 18 0 2 com- 
petitive isotope-labeling technique method of Payne et 
al [36]; however, on the basis of a steady-state analysis 
of their kinetic scheme, they argue that the rate of R 56 
directly depends upon the ratio of isotopically labeled 
C0 2 produced in their reactor. The large discrepancy 
between these studies has not been thoroughly explored 
and remains a bit of a mystery. Indeed, all three room 
temperature determinations may have been affected by 
unforeseen problems with handling the reactants or by 
unforeseen subtleties and complications in the overall 
reaction systems used. A glimmer of this is perhaps 



Figure 6 Low-pressure limits rate coefficient data of R 31 . 
Symbols are as follows: plus, Nettleton and Stirling [45]; 
open square, Westenberg and deHaas [16]; open triangle up, 
Astholz et al. [18,19]; open triangle right, Atkinson and Pitts 
[17]; filled triangle left, Merryman and Levy [44]; open dia- 
mond, values from the rate coefficient expression of Smith et 
al. [20]; star, Naidoo et al. [22]; open circle, this study. Lines 
are (in cm 6 mol -2 s -1 ) long dashed, Troe [21], & 3 i,o/[Ar] = 
4.0 x 10 28 T - 40 exp(— 2650 K JT), T = 300-2500 K; 
dot-dashed, Naidoo et al. [22], & 3 i,o/[Ar] = 2.4 x 10 27 
T~ 3 6 exp(— 2610 K/T), T = 220-2500 K; solid, this study, 
fc 3 i,o/[Ar] = 2.9 x 10 35 T - 6 0 exp(-4780 K/T) ± 6.1 x 
10 24 y— 3.0 exp (_i980 K/T), T = 300-2500 K. The aver- 
age deviation of the high temperature experimental data of 
Astholz etal. (T > 1700 K) [18,19] from our fit is ca. ±30%. 

the negative reaction rate observed by Graham et al. at 
their higher S0 2 concentration. To date, the generally 
accepted community value for R 56 is that of Graham et 
al. recommended by the IUPAC Subcommittee on Gas 
Kinetic Data Evaluation for Atmospheric Chemistry 
[35]. 

Recently, there have been two theoretical investiga- 
tions of R 56 going to products relevant to combustion 
conditions. Hindiyarti et al. [39] performed ab initio 
energy and structure calculations coupled to transition 
state theory (TST) to obtain a rate coefficient expres- 
sion for k- 56 (S0 2 ± H0 2 <— SO 3 ± OH) given as 
4.8 X 10 4 T 2 - 46 exp(— 13720 K JT) cm 3 mol " 1 s ” 1 in 
the temperature range of 800-2000 K. Wang and Hou 
[40] investigated the complicated reaction pathways 
leading from S0 2 ± H0 2 to the various products and 
posited that the effectively barrier-free HOSO ± 0 2 
(R_ 47 ) channel is the dominant reaction at “normal” 
temperature. Their argument can be presented in a sim- 
ple form as the competitive play between reaction dy- 
namics and reaction energetics: R_so (HS0 2 ± 0 2 «— 
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SO 2 + HO 2 ) is both kinetically and thermochemically 
disfavored (high barrier height and high endothermic- 
ity), R 56 is kinetically disfavored but thermochemically 
favored (a moderate barrier but strong exothermicity), 
and R _47 is kinetically favored but thermochemically 
disfavored (no effective barrier but moderate endother- 
micity). Based upon this argument, only the rate coef- 
ficients for R _47 (and stabilization of the intermediate 
complex, important near room temperature) were cal- 
culated using TST and a master equation analysis. The 
results pass through the middle of the room temperature 
experimental data and are claimed to support the IU- 
PAC evaluation. It should be noted that the difference 
between the calculated endothermicity of R _47 and the 
calculated barrier height for R 56 is only 35 kJ mol -1 
for their RCCSD(T)-cf calculation. Furthermore, while 
they do a good job of calculating the endothermicity of 
R —50 and R 56 they may be low for R _47 (see their Table 
1). As both R _47 and R 56 are competitive channels that 
pass through the same intermediate complex, increas- 
ing the endothermicity of R _47 increases the relative 
flux through the R 56 channel making it a significant 
reaction pathway. Even so, accepting Wang and Hou’s 
basic argument we would still expect no net flux, and 
thus no significant effect upon the course of reaction, 
for the R —47 channel given the excess O 2 in our system 
and the barrier-free reverse reaction back to reactants 
for that channel. This was borne out in our simulations 
where we found that contribution to SO 2 oxidation by 
either R 47 (R- 47 ) or R 50 (R- 50 ) is insignificant un- 
der our experimental conditions ( T = 970-1150 K). 
Rather, R 56 plays an increasingly important role as the 
temperature is lowered and the SO 2 concentration is 
raised. For instance, using our highest SO 2 mixture, 
S0 2 /H 2 /0 2 /Ar = 3.0/0.5/10.0/86.5, at T = 1033 K 
and P = 1.45 atm, ca. 32% of S0 3 was formed by R 56 , 
whereas ca. 66 % of SO 3 was formed by R 31 at times 
between 0.5 and 2.0 ms; while for our lowest SO 2 mix- 
ture S0 2 /H 2 /0 2 /Ar = 1.0/0.5/10.0/88.5 at T = 996 K 
and P = 1.31 atm, ca. 16% of SO 3 was formed by R 56 , 
whereas ca. 82% of SO 3 was formed by R 31 at the same 
time interval. The optimized kse values are plotted in 
Fig. 5 together with those of Hindiyarti et al. calcu- 
lated from K c , 56 , and the various room temperature 
measurements. Our data are well fitted as 

k 56 = 1.36 x 10 11 exp(— 3420 K/ T) cm 3 mol -1 s -1 

at T = 970-1150 K with estimated error limits of 
zb 30% and is shown as the dark solid line, which is 
extrapolated to room temperature (dashed line). It is 
interesting to note that the room temperature value of 
our expression is not far above the upper limit value 
of Graham et al. [37], the value recommended by the 


IUPAC Subcommittee on Gas Kinetic Data Evaluation 
for Atmospheric Chemistry [35]. As an enlightening 
exercise, we can find an equally good fit to our data 
by employing a gentle upwardly curving expression 
(of the basic form predicted by Wang and Hou [40]) 
that passes through the data of Payne et al [36]. This is 
shown as the light solid line in Fig. 5. Indeed, due to 
experimental scatter and the relatively narrow temper- 
ature range our data appear as essentially a single point 
on a graph like Fig. 5. Thus we can always find a good 
curve fit that passes through both our data and nearly 
any given room temperature value. Our data cannot 
help elucidate the rate of this reaction at lower temper- 
atures. As we have noted earlier, this channel becomes 
more important with decreasing temperature and so 
can play an important role in combustion exhaust. 

Our &3i ; eff values show, first, a very slight nega- 
tive temperature dependence at T = 970-1150 K and 
p ave = 16.2 pmol cm -3 and second, about 15-30% 
higher than those from & 3 i ? o expression of Naidoo et al. 
[22] at the same conditions. This is the indirect in- 
dication that we obtained the low-pressure limit rate 
coefficients fei 5 o). Thus, at each experimental condi- 
tion, k 3 i,o/[M] was obtained by scaling £ 31,0 with the 
corresponding collision-efficiency-corrected density, 
[M] = p x [1 + xo 2 (so 2 - 1) + *so 2 (£so 2 - 1)]. Here 
the collision efficiencies are £02 = 1.3, £so 2 = 9.0 and 

Mother species = 1-0, including £ A r- The & 3 i,o/[Ar] values 
obtained in this manner are listed in Table 3 in the 
Supporting Information and plotted in Fig. 6 . 

All possible uncertainties in our k 3 i,o/[Ar] values 
are estimated as follows: The presence of H atoms 
from sources other than the initiation reaction, R_n 
(HO 2 + H H 2 + O 2 ), and the uncertainty in [S 02 ]o 
are considered first. The effect of H-atom impurity 
was explored by injecting [H ] 0 = 0.1 ppb into the 
system — this [H]o is about five times larger than the 
kinetically produced [H]o in our experimental condi- 
tions. The ^310 values were increased by less than 20% 
and 3% with 1.0% SC^/^/CVAr mixture and 3.0% 
S 02 /H 2 / 02 /Ar mixture, respectively. The effect of the 
0.5% uncertainty in the nominal [S 02 ]o mole fractions 
was negligible. Here we placed a factor of 1.2 un- 
certainties in 10 owing to the possible presence of 
impurity as H atoms. 

The influence of the uncertainties from other reac- 
tion rate coefficients on the determination of i,o/[Ar] 
is examined next. As seen in the sensitivity figure, R 9 
(H + 0 2 + M -> H0 2 + M) exhibits the largest sensi- 
tivity. For the rate coefficients of R 9 we used k 9 ? o/[Ar] = 
5.55 x 10 18 r ~ 115 cm 6 mol -2 s -1 of Hwang et al. 
[24] with uncertainty limits of 30%. The uncertainties 
in k 3 i 5 o/[Ar] due to the coupling between R 31 and R 9 
were estimated by varying the rate coefficients of R 9 of 
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Hwang et al. by the given uncertainty limits and then 
observing the propagated effects on ^3i,o/[Ar]. With 
30% reduction of &9,o/[Ar], the resulting changes in 
&3i,o/[Ar] are 90%, 50%, and 40% for 1.0%, 2.0%, and 
3.0% S0 2 /H2/0 2 /Ar mixtures, respectively. But, for 
1.0% S 0 2 mixtures, the changes in reaction profiles 
were so great that reasonable matches to the character- 
istic times with any combination of k3i,o/[Ar] and k 56 
could not be achieved satisfactorily. We put a factor of 
1.6 uncertainties on &3i,o/[Ar] due to the coupling of 
R31 with R9. In the presence of S 0 2 , H atoms formed 
in the H 2 / 0 2 system are also consumed by R29 (S 0 2 ± 
H (±M) -> HOSO (±M)). The rate coefficient expres- 
sions of R 2 9 (&29,o/[M] and k 2 9 ,oo) are those of Blitz 
et al. [43]. We used our conservative error limit es- 
timate of ±30% for &29- The ±30% variation of k 2 9 
resulted in less than 20% change of k3i,o/[Ar]. So we 
set a factor of 1 .2 uncertainties on &3i,o/[Ar]. The effect 
of R14 (H 0 2 ± OH H 2 0 ± 0 2 ) on k 3 i, 0 /[Ar] via 
R56 (S 0 2 ± H 0 2 — > SO3 ± OH) was ca. 6% when ku 
values were reduced by ca. 35%. Here we set a factor 
of 1.1 uncertainties on k3i,o/[Ar] (see details below). 

The differences between the experimental and com- 
puted characteristic times for individual runs are on 
average less than 2% for 150 and 175 and 4% for r 2 5 
and T max . Taking into account all of the uncertainties 
discussed above plus additional uncertainties for Po , 
shock speed, signal-to-noise ratio, etc., we estimated 
maximum uncertainty limits, U = ^EI/ 2 = 3.0 for 
individual runs. 

A good fit to the data of this study and the previous 
experiments [16-19,22,44,45] is found using a combi- 
nation of two non- Arrhenius-type expressions given as 

^3i,o/[Ar] = 2.9 x 10 35 7’ _6 '° exp(— 4780 KJT) 

+ 6.1 x 10 24 7’ -30 exp(— 1980 KJT) cm 6 
x mol -2 s -1 

at T = 300-2500 K. As seen in Fig. 6, our data show 
the largest scatter. The rate coefficient expression 
shown above represents all the data within 26% error 
limits. Although we were unable to perform this study 
using N 2 as a collider gas, those modeling this reaction 
in Air might want to use the commonly referenced 
collider efficiency ratio of 0.77 (Ar/N 2 ) recommended 
by Atkinson et al. [35]. 

The rate coefficients from our expression are com- 
pared with those of previous studies. At T = 1050 K, 
the rate coefficient (in unit of cm 6 mol -2 s -1 ) from 
the above equation (3.1 x 10 15 ) is about a factor of 7 
and 3 smaller than the extrapolated value of Astholz 
et al. (2.2 x 10 16 ) [18,19] and Smith et al. (9.0 x 10 15 ) 
[20], respectively. But it is about 15% larger than the 


fit of Troe (2.7 x 10 15 ) [21] and RRKM fit of Naidoo 
et al. (2.7 x 10 15 ) [22]. There is a competition for S0 2 
between O and H0 2 . If H0 2 loss due to Ri 4 (H0 2 ± 
OH — >► H 2 0 ± 0 2 ) were increased, then there would 
be more S0 2 available and the overall reaction would 
be delayed even more. This would then require a de- 
crease in k3i ; o/[Ar] to account for the observed char- 
acteristic times, particularly at both low [S0 2 ]o and 
low-temperature conditions. The reverse argument is 
also true, and there is a range of ku values to choose 
among (see below). Using the lowest reported ku ex- 
pression, Kappel et al. [46], instead of Srinivasan et al. 
[25] used in this study, only increases the k3i,o/[Ar] 
values by 6% at our most sensitive conditions. 

The situation for R14 is certainly confusing at this 
time as can be seen in the rather busy nature of 
Fig. 7. The reaction has been studied quite extensively, 
both directly and indirectly, both experimentally and 
computationally, by numerous techniques covering the 
temperature range relevant from atmospheric chem- 
istry to flames. There is general agreement that the 
rate coefficient values of Keyser [47] at low tempera- 
tures and Goodings and Hayhurst [48] at high temper- 
atures are firmly rooted; for the midtemperature range 



Figure 7 Rate coefficients of R 14 (H0 2 ± OH -> H 2 0 + 
0 2 ). Symbols are as follows: plus, Cox et al. [56]; open trian- 
gle down, Sridharan et al. [57]; star, Goodings and Hayhurst 
[48]; open diamond, Keyser [47]; open square, Hippier et al. 
[50]; filled circle, Kappel et al. at P = 4 bar [46]; filled trian- 
gle up, Srinivasan et al. [25]. Lines are (in cm 3 mol -1 s -1 ): 
solid connecting two open diamonds, Hippier and Troe [49]; 
dashed, Gonzalez et al. [54]; solid, Kappel et al. [46]; solid 
connecting two open triangles up, Srinivasan et al. [25]; dot- 
dashed, Sivaramakrishnan et al. [52]; light solid, Rasmussen 
et al. [51]. [Color figure can be viewed in the online issue, 
which is available at www.interscience.wiley.com.] 
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of interest in this study, the values are still in flux. 
Troe and coworkers [46,49,50] have published a se- 
ries of increasingly refined shock tube measurements 
of R14 using various optical diagnostics to measure 
OH, HO 2 , and H 2 O 2 temporal profiles. As can be seen 
in Fig. 7, their earlier data (930-1680 K) [50] had 
a deep and narrow “V” temperature dependence cen- 
tered around 1250 K (open square). Their later data 
[46] had a similar but less pronounced feature cen- 
tered around 1000 K (filled circle). Nevertheless, they 
reported rate coefficient expressions that were either 
a straight line, ku = (2 ± 1 ) x 10 13 cm 3 mol -1 
s — 1 at T — 1000-1250 K [49] (solid line connect- 
ing two open diamonds), or an imperceptibly curved 
line, k u = 1.0 x 10 13 + 5.8 x 10 13 exp(-2000 K JT) 
cm 3 mol - 1 s - 1 at T = 950-1250 K [46] (dark solid 
line). More recently, Srinivasan et al. [25] performed 
a direct determination of £14 values via a joint opti- 
mization of the reactions OH + NO 2 -> HO 2 + NO 
and H0 2 + OH H 2 0 + 0 2 (R14) reactions behind 
reflected shocks using an OH diagnostic. Their recom- 
mended rate coefficient expression is ku = (3.0 zb 1.8) 
x 10 13 cm 3 mol -1 s _1 at T = 1200-1700 K (solid 
line connecting two open triangles up) — this ku ex- 
pression was used throughout this study. This expres- 
sion is 25% lower than the best fit to their data, but 
they concluded that it is the best available experimen- 
tal value for combustion-modeling applications based 
upon their consideration of other high-temperature 
data. Rasmussen et al. [51] performed a critical review 
of all the available data for this reaction as part of the 
analysis of their high-pressure laminar flow measure- 
ments of the C 0 /H 2 0 / 02 /N 0 x reaction system. In their 
temperature range, 600-900 K, there are no experimen- 
tal measurements of R 14 . They proposed an expression 
that is the sum of three three-parameter rate coefficient 
expressions that does an excellent job of fitting the 
high-temperature data of Goodings and Hayhurst and 
Kappel et al., including the sharp dip, as well as the 
low-temperature data of Keyser. It also passes through 
the data of Srinivasan et al. (filled triangle up). Their 
expression is a factor of 3 lower than ours on average 
in the range of mutual overlap. Sivaramakrishnan et al. 
[52] derived an expression for ku indirectly as a result 
of their high pressure single pulse shock tube study of 
the oxidation of diluted CO mixtures doped with small 
quantities of H 2 . Their expression is given as the sum of 
five three-parameter rate coefficient expressions over 
the temperature range 298-1600 K that produces a nar- 
row dip at T = 1400 K (dot-dashed line). They argue 
that this expression is needed to match their CO and 
CO 2 concentration data at their highest temperatures, 
1400-1500 K, and their highest pressures, 256 and 450 
atm. The dip in their expression at 1400 K is shifted 


upward in temperature with respect to Hippier et al. 
(1250 K) and Kappel et al. (1000 K); it is also neither as 
deep nor as sharp. The Sivaramakrishnan et al. and our 
ku expressions agree within 10 % over the temperature 
range of mutual overlap. Their expression decreases 
with temperature in this range, and they are mostly 
above us with the expressions crossing ca. 1075 K (the 
upper third of our temperature range). The most recent 
theoretical study of R 14 has been performed by Gonza- 
lez et al. [53,54] (see references therein for a thorough 
review of previous theoretical work), where they deter- 
mined ab initio geometries at the HF and MP2 levels 
and energies at the MP4 level. They are careful to de- 
tail the problems with calculating the surfaces at this 
level of theory. For the singlet reaction surface, they 
determined that this pathway would be unimportant 
near room temperature, their area of focus, due to the 
63.6 kJ mol -1 barrier heights. For the triplet reaction 
surface, they do predict that reaction can occur at room 
temperature. However, they cautioned that their calcu- 
lations were not sufficiently adequate to quantitatively 
determine the coupling between the nonplanar surfaces 
at the MP2 level. Given that caveat, they moved for- 
ward to predict a rate coefficient expression using vi- 
brational/rotational adiabatic theory. Their calculated 
temperature dependence of ku values was then com- 
bined with the room temperature experimental data to 
produce the best estimate of the ku expression given as 
£14 = 4.3 X 10 13 (77300)-°- 21 exp(57 K IT) cm 3 mol " 1 
s-' at T = 298-1 100 K (dashed line). 

The direct experimental evidence has produced 
three effectively straight line rate coefficient expres- 
sions, each within the error bars of the others, all lying 
within the range (2-3) x 10 13 cm 3 mol -1 s -1 and over- 
lapping around 1250 K. It cannot be overemphasized 
that none of them shows any evidence of the narrow 
“V” behavior. The indirect experimental evidence has 
produced a very complex expression but one that is 
within 10% of the Srinivasan et al. expression used 
in this study around 1000 K and well within the error 
bars of all three direct studies. The available theoret- 
ical evidence, while carefully qualified by its authors, 
gives a very mildly curved rate coefficient expression 
that passes through the error bars of the direct experi- 
ments. Combining the published data given in the three 
studies with the low-temperature data of Keyser and 
the high-temperature data of Goodings and Hayhurst 
would strongly suggest the same sort of broad and shal- 
low bowl- shaped appearance that is the signature of the 
HO 2 self-reaction. Clearly, the advances made by the 
theoretical community during the intervening decades 
since the most recent computational study would help 
to elucidate the nature of this important reaction and 
guide future experiments. 
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Figure 8 (a) Conversion yield of S(IV) to S(VI), calculated at a constant density condition (p = 16.0 pmol cm -3 ) using the 

Table 2 mechanism as functions of temperature and reaction (residence) time. Symbols are for reaction time: triangle, 1 ms; 
square, 2 ms; circle, 5 ms. (b) Conversion yield calculated at a constant pressure condition (P = 1 atm). Symbols are as in 
panel a. 


Among our mixtures, 1.0% S0 2 gave the highest 
conversion yield of S(IV) to S(VI), i.e., S0 2 to SO 3 (the 
contribution to S(VI) from H0S0 2 and H 2 SC >4 being 
negligible at our experimental conditions). Thus we 
computed the maximum conversion, [S 03 ] max /[S 0 2 ]o, 
in the S0 2 /H 2 /0 2 /Ar = 1.0/0.5/10.0/88.5 mixture us- 
ing the Table 2 mechanism as functions of temperature 
and reaction time (residence time) at a constant density 
(p = 16.0 pmol cm -3 ) and a constant pressure (P = 
1.0 atm) flow conditions. As seen in Figs. 8 a and 8 b, 
for either flow condition, maximum conversion occurs 
at temperatures around 1200 K and slightly increases 
with reaction time. The maximum conversion yield is 
ca. 3.5% and 2.5% for the constant density and con- 
stant pressure flow condition, respectively. Flux anal- 
yses show that SO 3 is mainly produced by R 31 (S0 2 
+ O + M S0 3 + M). At low temperatures, there 
is some contribution from R 56 (S0 2 + H0 2 — > SO 3 + 
OH). The major consumption channels are R 34 (SO 3 
+ H — > S0 2 + OH), R 35 (SO 3 + O — > S0 2 + 0 2 ) and 
R_ 5 2 (HOS0 2 + 0 2 SO 3 + H0 2 ). It is noteworthy 
that the maximum conversion occurs at temperatures 
above that of the maximum ^310 values. Undoubtedly, 
this is because that the conversion yield is dependent 
upon not only the rate coefficient values of £ 31,0 and 
&56 but also the availability (concentrations) of H, O, 
and H0 2 in the system. 

The unusual behavior in the rate coefficients of the 
spin-forbidden title reaction, R 31 (S0 2 (X 1 Ai) + 0( 3 P) 


+ M S0 3 (X 1 A / 1 ) + M) has been rationalized by 
Astholz et al. [19] in terms of the potential barrier aris- 
ing from triplet-singlet crossing of the SO 3 complex 
formed (S0 3 ( 3 A / ) SC^A')) and decrease of the 

collision efficiency (/3 C ) with temperature in stabilizing 
the collision complex, i.e., the rate coefficients are con- 
trolled by the competition between the barrier crossing 
and the stabilization of the complex. Astholz et al. used 
Vo = 13.8 it 4 kJ mol -1 and a factor of 4 decrease of 
from 300 K (/ % = 0.4) to 2000 K ( = 0.1), whereas 
Naidoo et al. [22] used 15.9 kJ mol -1 in their calcu- 
lations. However, for understanding of this reaction in 
more detail, high-level quantum calculations are bene- 
ficial to elucidate the nature of the triplet SO 3 potential 
curve — it may be repulsive around the crossing — and 
the coupling between triplet-singlet states. Recogniz- 
ing the triplet-singlet crossing actually occurs along 
the seam at the lowest crossing point, it is highly de- 
sirable to calculate the potential energy surface (PES) 
around the crossing point. If the seam crossing prob- 
ability is energy dependent, a new temperature effect 
could be introduced. 

We now consider the possibility of SO 3 formation 
by a molecular reaction channel involving the relatively 
abundant amounts of S0 2 and 0 2 available in the hot 
oxidizing environment of the high-pressure section of 
an aircraft engine. The most efficient version of the 
reaction would be the direct oxidation of S0 2 to SO 3 
in the form of R 58 (S0 2 + S0 2 + 0 2 — >- SO 3 + SO 3 ). 
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Such a reaction may be viewed as proceeding, similar 
to S20^ reactions [ 55 ], via SO2 + O2 ^ SO4 + 
SO2 — O3SOSO2 — SO3 SO3 or SO2 + SO2 
0 2 0 2 S 00 S 0 2 -> S 0 3 + S 0 3 . Compared to S 2 Og 

fragmentation reactions, either reaction shown above 
does not have an available donor electron for easy frag- 
mentation to SO3. Also a large number of rearrange- 
ments for intermediate formation would be required, 
so that high reaction rate coefficient values are not ex- 
pected. Nevertheless, we can assign a maximum value 
to the rate coefficient of R58 at our test conditions based 
upon the observed flatness of the absorption profile 
preceding the exponential rise due to chain reaction. 
The proposed molecular reaction would start at time 
zero and be evidenced by a steady decrease in the ini- 
tial absorption, entirely due to SO2, as expected for a 
mixture containing both the reactants. The k 53 value 
was chosen such that visible change in the initial ab- 
sorption signal, greater than our signal-to-noise ratio 
was noticed. For our conditions, £53 = 3.0 x 10 12 cm 6 
mol -2 s -1 was chosen (the reverse rate coefficient was 
calculated via K c ^). First, we chose the highest SO2 
mixture run at the lowest temperature and second, be- 
cause the influence to the conversion yield from S(IV) 
to S(VI) by R58 would increase with the increase of 
[S 0 2 ]o and decrease with T, we selected a test condi- 
tion of T = 1050 K and p = 16.0 pmol cm -3 in 3 . 0 % 
SO2 -containing mixture. As seen in Fig. 9 , the charac- 



Figure 9 Simulated absorption profiles using the Table 2 
mechanism with and without R 58 (SO2 + SO2 + 02^ 
SO3 + SO3) for the condition of S02/H2/02/Ar = 
3.0/0.5/10.0/86.5, T = 1050 K and p = 16.0 ^mol cm" 3 . 
Lines are solid, without R 58 ; dot-dashed, with R 58 and 
&58 = 3.0 x 10 12 cm 6 mol -1 s -1 ; dashed, with R 58 (£ 5 8 = 
3.0 x 10 12 cm 6 mol -1 s -1 ) and £31,0 (SO2 + O + M -> 
S 0 3 +M) x 1.26. 


teristic times (125, 150, 175, and r max ) were shortened by 
ca. 13 % and the maximum normalized growth constant 
of the absorption signal ([dA(t)/dt/A m3LX ] m3LX ) — shape 
of the absorption signal — was increased by ca. 20%. 
With the increase of only the ^3 10 value by 26 % (still 
within our stated upper error limit), satisfactory agree- 
ment was achieved in both the characteristic times and 
the maximum normalized growth constant. At high 
temperatures for 3 %, 2 %, and 1 % SO2 -containing mix- 
tures ( 1150 , 1150 , and 1040 K, respectively), however, 
the simulated absorption profiles with inclusion of R58 
showed no alteration of the reaction characteristics. 
The changes in the low temperature {t;} values are 
due to reduction of reaction flux through R31 (SO2 + 
O + M S 0 3 + M) by depletion of available [S 0 2 ] 
by R 58 (S 0 2 + S 0 2 + 0 2 -* S 0 3 + S 0 3 ). 

With the inclusion of R58, as expected, we observed 
some increase of the conversion yields at all temper- 
atures for the 3 % and 2 % SO2 -containing mixtures. 
For example, with the rate coefficient and test condi- 
tions given above the conversion yield was increased 
from 1 . 7 % to 2 . 5 % at 5 ms where the maximum con- 
version occurs. We also observed the increase of the 
conversion yield at all temperatures for the 3 % and 2 % 
SO2 -containing mixtures. 

It is apparent that R58 affects the characteristic times 
and conversion yields for high [S02]o-containing mix- 
tures. But the extent of the effect is not large: As shown 
above, agreement in the characteristic times with and 
without R58 could be achieved by increasing ky\^ value 
within its error limit. And the absolute values of the 
conversion yield with R58 are still small. Lack of exper- 
imental or theoretical evidence for R58 and the reasons 
given above permitted us to exclude R58 in the Table 2 
mechanism. 


CONCLUSIONS 

Thermochemical and structural data for HOSO, HSO2, 
and HOSO2 were critically reviewed, and the ther- 
modynamic data for these species were estimated and 
NASA-type polynomial coefficients were established 
(Table 1 in the Supporting Information). 

The title reaction was studied in the perturbed sys- 
tems where the well- studied F^/C^/Ar system was per- 
turbed with small amounts of SO2. A comprehensive 
reaction mechanism for SO2/H2/O2MX system was as- 
sembled (Table 2 in the Supporting Information). Op- 
timum S02/H2/02/Ar mixture compositions and ex- 
perimental conditions were obtained via a series of 
sensitivity analyses. Shock tube experiments were per- 
formed using an OH laser absorption spectroscopic 
technique. The measured absorption coefficients of 
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50 2 at \q = 310.032 nm (air) and T = 300-1150 
K are 

e c (S0 2 ) = 1.3 x 10 5 exp(— 180K / 7)cm 2 mol -1 

The experimental results were analyzed via computer 
simulations, in that the rate coefficients of R 31 (S0 2 + 
O T- M — > SO 3 + M) and R 56 (S0 2 -T H0 2 — > SO 3 -F 
OH) were simultaneously optimized against the exper- 
imental reaction characteristics (characteristic times). 
The rate coefficient expressions are 

k 56 = 1.36 x 10 11 exp(— 3420K/T)cm 3 mol -1 s -1 

at T = 970-1150 K with estimated error limits of ± 
30%, and 

* 3 i,o/[Ar] = 2.9 x 10 35 T~ 60 exp(— 4780 K/T) 

+ 6.1 x 10 24 T “ 3 0 exp(-1980K/r) 
x cm 6 mol -2 s -1 

at T = 300-2500 K and p ave = 16.2 pmol cm -3 with 
error limits of 26%. 

The calculated maximum S(IV) to S(VI) conver- 
sion yield, ca. 3.5% and 2.5% for the constant density 
and constant pressure flow condition, respectively, oc- 
curred at rather higher temperatures (M200 K) than 
lower where the maximum &310 value is located (~800 
K). We found that the conversion yield is controlled by 

50 3 production flux mainly through R 31 (S0 2 + O + 

M — > SO 3 + M) and, to the lesser extent, through R 56 
(S0 2 + H0 2 S0 3 + OH) and consumption flux 

through the reactions with H (R 34 , SO 3 + H -> S0 2 + 
OH), O (R 35 , SO 3 + O — > S0 2 + 0 2 ) and H0 2 (R_ 5 2 , 
HOS0 2 + 0 2 S0 3 + H0 2 ) as well. 


All experiments and computations in this study were per- 
formed at the NASA Glenn Research Center in Brook Park, 
Ohio. 
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